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Abstract 



We construct the quasi- stationary (QS) probability distribution for the 
Domany-Kinzel stochastic cellular automaton (DKCA), a discrete-time 
Markov process with an absorbing state. QS distributions are derived at 
both the one- and two-site levels. We characterize the distribuitions by their 
mean, and various moment ratios, and analyze the lifetime of the QS state, 
and the relaxation time to attain this state. Of particular interest are the 
scaling properties of the QS state along the critical line separating the active 
and absorbing phases. These exhibit a high degree of similarity to the contact 
process and the Malthus-Verhulst process (the closest continuous-time analogs 
of the DKCA), which extends to the scaling form of the QS distribution. 
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I. INTRODUCTION 



The Domany-Kinzel stochastic cellular automaton (DKCA) Pj] is a Markov process that 
exhibits a phase transition from an active state to an absorbing one. Stochastic processes 
with an absorbing state arise frequently in statistical physics @, and are currently of great 
interest in connection with self-organized criticality and nonequilibrium critical phenom- 
ena ||||. Many studies of the DKCA and related probabilistic cellular automata (PCA) 
have been published, using deterministic mean- field equations Monte Carlo simula- 

tions P H13I1 and renormalization group (RG) analyses []14|-|T^ . 

While the mean-field (MF) description of the DKCA admits (for appropriate parame- 
ter values) an active stationary state, for finite system sizes the model always ends up in 
the absorbing state, due to fluctuations. MF theories ignore such fluctuations, and so are 
incapable of treating finite systems. But, since simulations and other numerical methods 
typically study finite systems, it is of interest to develop approximate theoretical descriptions 
that account for finite system size. A natural way to study finite systems with an absorbing 
state is via the quasi- stationary distribution, which, when it exists, describes the asymptotic 
properties conditioned on survival ]TT|-|TU||. Recently, mean-field-like methods were devel- 



oped for studying the quasi-stationary state of finite systems with an absorbing state [20 



The quasi-stationary properties converge to the true stationary properties in the infinite- 
size limit. (Indeed, this provides the rationale for studying the "stationary" behavior of 
absorbing-state models in simulations, which of necessity treat finite systems.) In Ref. [GO 



quasi- stationary distributions for various continuous-time Markov processes are constructed, 
in particular, for the contact process (CP) and the closely related Malthus-Verhulst process 
(MVP). In the case of the CP, both one- and two-site approximations are derived. In this 
work, we extend the analysis to discrete-time processes, using the DKCA as an interesting 
example, closely related to the CP. 

The remainder of this paper is structured as follows. In Sec. II we review the one- and 
two-site MF approximations for the DKCA. In Sec. Ill, we construct the quasi-stationary 
(QS) probability distribution at the site level. The QS distribution at the pair level is 
discussed in section IV, while Sec. V presents our conclusions. 



II. SITE AND PAIR MF APPROXIMATIONS 

In this section we review the definition of the DKCA and its mean-field description at 
the one- and two-site levels The DKCA is a discrete-time Markov process (all sites 

are updated simultaneously), whose configuration is given by a set of stochastic variables 
{o~i} = or 1), defined at sites i and times t = 0, 1, 2,..., such that t + i is even. Let a 
represent the configuration at time t, and -Pt(cr) the probability distribution in configuration 
space. The evolution of the latter is governed by 

P t+1 (a) = J>Ha')P,(a') (1) 

cr' 

where u(a\a') denotes the probability of the transition a' — > a, and enjoys the properties 
oo(a\a') > and Y^a u(a\a') = 1. The transition probability for the DKCA is a product of 
factors associated with each site: 
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u{a\a') =Hwi(a i \a r ) , (2) 

8=1 

where w^a^a') is the conditional probability for site % to be in state at time t + 1, given 
configuration a' at time t. The probabilities Wi(<Ti\a') are translation- invariant and in fact 
depend only on the variables o\_ x and cr' i+1 at the previous time step: 

Wi(oi\(r') = w D * : (cr i |o--_ :l ,er- +1 ) . (3) 

The above relations, with the transition probabilities given in Table I, define the DKCA. 
Noting that the transition (00) — > (1) is prohibited, we see that the configuration <7j = 0, Vi 
is absorbing. 

Of interest are the n-site marginal probabilities. The evolution of the one-site distribution 
Pt(o-i) = Eajj^iPtio-), is given by, 

Pt+i(o-i) = X X w DK{(Ti\o-' i _ 1 ,a' i+1 )P t (a' i _ 1 ,a' i+1 ) , (4) 

where P^(a'^_ 1 , is the marginal distribution for a pair of nearest-neighbor sites. The 
evolution of the latter is coupled to the three-site probability, so: 

P t +i(o-i-i, a i+ i) = X XX ^)«W^+ik', ^+2)^(^-2) ^+2) • (5) 

Evidently we have an infinite hierarchy of equations. In the n-site approximation the hier- 
archy is truncated by estimating the (n+l)-site probabilitites on the basis of those for n 
sites. 

The simplest case is the one-site approximation H||], in which P t {a[_ 1: cr,- +1 ) is factored 
into a product of one-site probabilities. This yields the recurrence relation, 

Pt+i = Pt [2pi - (2p 1 -p 2 )p t ] , (6) 

where pt = Pf(l) is the density of active sites (the order parameter for the DKCA). Eq. 
(P) admits two stationary solutions, corresponding to the possible DKCA phases: absorbing 
(p = 0), and active, in which, for p\ > 1/2, 

P=^. (7) 
2jJi - P2 

Thus the critical line at the site level is p\ c = 1/2. 

In the pair approximation |J the three-site probability is written in terms of the two-site 
quantity, using the conditional probability: 

Pt(0~i) 

(The one-site probabilities are given by Pt(<7j) = J2cr i+2 Pti.°~ii °~i+2))- Letting z t = P t (l, 1), 
and noting that Pt(l, 0) = k t = pt — Zt, we have the relations, 



3 



k 2 

z t +i = — [P2Z t + Pik t } 2 +pl ' 



1 

Pt 



I-Pt 



and 



kt+i = —{p2Zt+pik t ){q2Zt+qih) + , Plkt 

Pt 1 - Pt 

where qi = l—pi, while for v t = P t (0, 0) we have: 



(qxkt + u t ) 



(9) 



(10) 



v t +i = — [q2Zt+qik t f + z ~ — + • (11) 
p< i - pt 

In the active stationary state, these relations imply, 

l-2pi n9 , 

« = 7T- P » 12 

P2 - 2pi 

which leads to the stationary active-site density, 

= P2 (p 1 -l) 2 + p 1 (3 Pl -2) 
P (2 Pl -l)(2 Pl -p 2 ) ' 1 j 

In this approximation, the critical line in the (pi,P2) plane is: 

Pi(2-3pi) 

P2 = -p. vj- • 14 

(i -pi) 1 

The phase diagram for the DKCA in the one- and two-site approximations is compared with 
simulation results 21] in Fig. 1. 



III. QUASI-STATIONARY PROBABILITY DISTRIBUTIONS 

A. Method 

Since the approximations discussed in the previous section effectively consider the L — > oo 
limit, the densities (p, z, etc.), are in fact deterministic variables. Our goal in this paper 
is to construct reduced stochastic descriptions of a finite system, in a manner analogous to 
that employed in deriving n-site approximations, and to determine the associated quasi- 
stationary properties. In this section we study the problem at the one-site level. Consider 
the DKCA on a ring of L sites. At the one-site level, the state of the system is specified by 
N t , the number of active sites at time t. Let Pt{N) (N = 0, L) be the probability to have 
exactly iV active sites at time t. The probability vector p t = [p t (0),p t (l), ...,p t (L)], satisfies 

Pt (N)=J2W(N\N>)p t „ 1 (N>) , (15) 

N' 

where W is the transition matrix, with W(N\N') representing the probability to have iV 
active sites at time t, given N' at time t — 1. 
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At the one-site level, the state each site is treated as an independent event. Given N' 
active sites at time t, our best estimate for the probability x of a given site to be active at 
the next time step is (see Eq. (||): 

x = y [2pi - {2p 1 -p 2 )y\ , (16) 

where y = N'/L. Thus the transition probabilitites in the 1-site approximation are: 

w(Nm = x " ■ (17) 

(Here we suppose that all configurations with the same number of active sites are equally 
probable, since there is no reason to prefer one such configuration over another at this 
level of analysis.) The one-site distribution pt{N) is therefore a superposition of binomial 
distributions with means x = 0, 1, ...,L, the weight of a given distribution depending only 
on the mean population at the previous step. 

Since (for L finite) the probability distribution will always evolve to the absorbing state, 
p(N) = 5jv,o ; h is of interest to study the quasi-stationary distribution p(N), defined as 
follows. We suppose that as t — > oo, the probability distribution, conditioned on survival, 
attains a time-independent form. This means for long times 

Pt(N)=A t S N>0 + S t p(N), (18) 

where the only time dependence lies in A t and St- Since the QS distribution p(N) is condi- 
tioned on survival, p(0) = 0. Adopting the normalization 

E p( n ) = !> ( 19 ) 

N=l 

S t in Eq. (0) represents the survival probability and A t = I — St the probability to have 
fallen into the absorbing state. The QS hypothesis is verified numerically below. Evolving 
the distribution in Eq. ( |18~D to the next time step, we have 

L 

p t+ i{N) = A t 5 Nfi + S t J2 W(N\N')p(N') 

N'=0 

= A t+1 5 Nfi + S t+1 p(N), (20) 
which implies that St+i = a St, where 

L 

a = 1 - £ W{0\N)p(N) . (21) 

N=l 

Restricted to the states 1, ...,N, p is an eigenvector of matrix W with eigenvalue a. The 
lifetime r of the QS state is: 

r = ~ ■ (22) 
ma 
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One method for generating the QS distribution is via iteration of the evolution equation, 
Eq. fll5D , until the distribution q t (N) = p t (N)/S t (for N = 1,2,...,L), attains a time- 
independent form. We refer to this as the direct method. An alternative method |23|] is 
based on writing the evolution in the form: 

Ap t (N)=p t+1 (N)-p t (N) = -w(N)p t (N)+r t (N), (23) 

where r t (N) = J2N^NW(N\N')p t (N') and w(N) = En>?nW(N'\N). Inserting the nor- 
malized QS distribution p(N) in the r.h.s. of the above relation, we have 

(a-l)p(N) = -w(N)p(N) + r(N), (24) 

where r(N) = Y^N'W(N\N')p(N'). Noting that I — a = r(0), this may be written in the 
form 

p(/v » = ^fW- n* 1 - (25) 

This relation suggests the following iterative scheme: 

P'W = MN) + {l-a) w( ^ r(oy (26) 

where a is a parameter and r(N) is evaluated using the distribution p(N). At each it- 
eration the new distribution p' must be normalized. In this way, we can construct the 
quasi- stationary state from any initial initial distribution p(N) that is nonnegative and nor- 
malized. We call this the iterative scheme. As discussed in Ref. [f23j, good convergence is 
obtained for a ~ 0. 



B. Results 

We have constructed the QS distribution for the DKCA at the one-site level, using both 
the direct and iterative schemes. In Fig. 2 we show the time evolution of the probability 
distribution (conditioned on survival) at a point on the critical line. It is evident that the 
distribution reaches a quasi-stationary form after about 100 time steps. Fig. 3 shows the 
evolution of the mean population (N), the moment ratio m = (N 2 ) / (N) 2 , (both conditioned 
on survival), and the decay rate 7 = r t (0)/S t , to their quasi-stationary values, for the same 
parameters as in Fig. 2. (Note that 7 is the transition rate into the absorbing state.) 

Relaxation to the QS state appears to consist of two stages: an initial transient, which 
depends strongly on the initial condition, and a long-time, exponential approach to the final 
values. The mean population, for example, follows \(N) — (N)qs\ ~ e~ t l TR . (For the data in 
Fig. 3, t r ~ 9, while the lifetime of the QS state is about 16.) We find that m and 7 relax at 
the same rate as (N), and that this relaxation time is independent of the initial distribution. 
Thus the asymptotic relaxation to the QS state is governed by a relaxation time tr that 
appears to depend only on the parameters p\ and p<i and on the system size L. 

The distribution may be further characterized by its skewness, S, and kurtosis, K, defined 
through the relations: 
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and 



K 2 



K = 5, (28) 



where K n is the n 4/l cumulant of the distribution M. For the Gaussian distribution, both 
skewness and kurtosis are null (S — K — 0). The evolution of the skewness and kurtosis 
is also shown in Fig. 3. Fig. 4 shows the quasi-stationary distribution at several points in 
parameter space. We observe that in the frozen phase the distribution collapses to N = 1 
while in the active phase it is concentrated near N = L. 

Of particular interest are the QS scaling properties at the critical point. We have verified 
that (N) ~ L 1 / 2 in the critical QS state. The QS lifetime scales in the same manner. These 
system-size dependences were encountered previously for the CP and MVP in the one-site 
approximation ||20|| . The relaxation time tr also grows ~ L 1 / 2 at the critical point; we find 



tqs/tr ~ 2.67 for Pl =p 2 = l/2. 

In the active phase, however, tqs grows ~ exp [const, x (p\ — pi c )L], while tr varies only 
slightly with pi, p 2 , and L. This leads to a clear separation of time scales (tqs ^> Tr) for 
large systems. (For L = 100, pi = 0.6, and p 2 = 0.5, for example, tq S — 2000 while t r ~8.) 

In Fig. 5 we show the QS density versus p\ in the site approximation, for several system 
sizes, showing convergence to the deterministic mean-field prediction. Also shown is the 
moment ratio m versus p\ for the same system sizes. Data for L = 1000 - 10 5 (see Fig. 6) 
indicate that as L — > oo at the critical point, the moment ratio approaches the value 1.660, 



found for the CP in the one-site approximation, and for the Malthus-Verhulst process |20 

The moment ratio m appears to approach the same limiting value all along the critical 
line, for p 2 < 1. This suggests that the critical QS distribution has a scaling limit for large 
L, of the form 

P(N) * ^-P(N/(N)) , (29) 

where V is a scaling function. Fig. 7 compares (for p 1 = p 2 = 0.5), (N) p(N) as a function 
of N/ (N), for system sizes L= 10 3 , 10 4 , and 10 5 , as well as the exact scaling function for the 
CP and MVP found in Ref. |20[|. It is interesting to note that the QS distribution for the 



DKCA has the same scaling form as for the CP and the MVP, despite the fact that in the 
critical DKCA, p(N) takes its maximum value for N > 1. The position of this maximum, 
however, grows very slowly with L (roughly, ~ In L) , so that it does not alter the infinite-size 
limit. A possible explanation for a maximum away from iV = 1 is that, in the DKCA, there 
are transitions to the absorbing state (N = 0), from various values of N, not only for N = 1, 
as is the case in the CP. 

The inset of Fig. 7 shows (N) p(N) versus N/ (N) for p\ = 0.5, L = 2000, and various 
values of p 2 . The data collapse confirms the scaling hypothesis, except for p 2 — 1. The 
distinct behavior in the latter case is expected, since P2 = 1 corresponds to compact directed 
percolation, which has two absorbing states (all or all 1). (The situation is analogous to 
that found in the voter model, a continuous-time process with two absorbing states [p0|] .) 
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We find that any initial distribution not concentrated on iV = evolves to the QS 
distribution, which is independent of the initial condition. (Uniqueness of the QS distribution 
is to be expected in the case of the DKCA, which has only one absorbing and one active 
state. A nonunique QS distribution can be envisioned for a process in which the active state 
exhibits symmetry breaking.) 

The direct and iterative methods discussed above yield (as they must), the same QS 
distribution. Using a in the range -0.4 to 0, we find that the iterative method converges to 
the QS distribution slightly faster than the direct approach (it typically requires about 30% 
fewer steps). In contrast with the continuous-time case, in which the iterative method can 
be orders of magnitude faster than integration of the master equation , here the gain 

in efficiency is quite modest. This is not surprising, since the enormous gain in efficiency 
for continuous-time processes is associated with the small time step required to maintain 
numerical stability, in the usual direct integration schemes. In the present case, the direct 
method has an effective time step of unity. 



IV. PAIR APPROXIMATION 
A. Method 

In this section we construct the quasi-stationary probability distribution for the DKCA 
at the pair level. The system is described by two stochastic variables, the number of occupied 
sites N, and the number of doubly occupied nearest-neighbor (NN) pairs, Z. We consider 
a ring of L sites. [For convenience we introduce a different notation for the site variables, 
defining ipi = <Ji/ 2 for even t and <fi = <J(i+i)/2 for odd t. In this way the site index always 
takes the values 1,...,L at all times, and NN sites have state variables ifii and y^.+i-] 

To begin, we establish the allowed range of values for Z. Using '1' and '0' to represent 
occupied and vacant sites, respectively, we denote by K the number of (10) NN pairs. (By 
symmetry, the number of (01) pairs is also K.) K is not an independent variable, since each 
1 is followed by a or another 1, yielding N = Z + K. Similarly, the number of (00) pairs, 
V, is given by V = L — 2N + Z. The conditions K > and V > imply certain limits for 
Zona ring of L sites, listed in Table II. 

Next we construct the transition probabilities W(N, Z\N', Z'). Note that the presence 
of V (00) pairs at time t implies that there are at least this many vacant sites at time t+1; 
thus W = for N > 2N' — Z'. We proceed by analogy with the one-site approximation: 
given N' and Z', we first determine the pair densities z = Z/L, k, and v using Eqs. (|9])- 
(|TT|). (Here z, k, and v represent the densities at time t + 1, while the variables appearing 
on the r.h.s. of each equation are evaluated using p t = N'/L and z t = Z'/L.) We treat all 
configurations having the same N and Z as equally probable, and estimate the probability 
of any one such configuration as: 



Q(N,Z;p,z) 



z Zf,2K v V 



p N {l-p) L ~ N 



jt 8 -^, (30) 
L j= i Pj 



where k = p+z. The first factor (in square brackets) is the product of all pair probabilities, 
divided by the product of all site probabilitites. The second factor represents a correction, 



S 



needed for normalization of Q, which arises as follows. Suppose we construct the probability 
Q starting at site j, so that the first factor in the product is pj, i.e., the pair probability 
associated with sites j and j + The next factor will then be pj + x/sj + i, which represents 
the conditional probability of the variable at site j + 2, given the state of site j + 1, and so 
on. When we close the ring, adding the link between sites j — 1 and j, the final factor is 
P(cpj-ij(pj\(pj-x,(fj) = 1, not P(ipj_i,(pj)/[P(ipj_i)P((pj)]. So the first factor in Eq. ( |30|) 
has one pair factor, and two site factors, too many, and should be multiplied by Sj-xSj/pj-i. 
Since the position of the starting link is arbitrary, we take the mean of this correction over 
the ring. 

Note that the correction factor may be written so: 
j_ s j s j+i _ ]_ 

Lfri p 3 L 

In case the pair numbers take their expected values (Z = Lz, etc.), the correction factor is 
unity. 

The transition probability is the product of a configurational probability Q(N, Z\ p, z) 
and the number of configurations, T(N, Z; L), having exactly N active sites and Z active 
pairs on a ring of L sites: 

W{N, Z\N', Z') = T(N, Z- L)Q(N, Z; p, z) , (32) 

for (N, Z) ^ (N f , Z') and N < 2N'-Z'\ for N > 2N'-Z', W = 0; if (N, Z) = (N 1 , Z'), 

W(N, Z\N', Z') = l- r(JV, Z- L)Q(N, Z- p, z) , (33) 

N,Z 

where (*) denotes the exclusion of the single term N = N', Z = Z'. An expression for the 
combinatorial factor T(N, Z; L) is derived in the Appendix. (Note that T = for values of 
A^ and Z outside the permitted range given in Table II.) The evolution of the probability 
distribution follows: 

p t+1 (N, Z)= ]T W(N, Z\N', Z')p t (N', Z') . (34) 

N',Z'\2N'-Z'>N 

B. Results 

We constructed the QS distribution at the pair level for the DKCA, for systems of up to 
200 sites, focusing on the behavior in the vicinity of the critical line. (The computation is 
considerably more demanding of memory and cpu time than is the one-site approximation; 
the chief limitation is the evaluation of the coefficients T.) 

Fig. 8 shows the QS order parameter versus p\ in the pair approximation, for several 
system sizes. We also plot the moment ratio m, showing a series of crossings whose location 
approaches the critical point as L — ► oo. The behavior of m versus L, at criticality, is shown 
in Fig 6. In Fig. 9 we show the QS distribution, p(N,Z), at criticality (j> 2 = 0.5,pi = 
0.6306), for L = 100. The marginal distribution p(N) is similar to that found in the site 
approximation. The behavior of the mean population, moment ratio, decay rate, skewness 
and kurtosis, as functions of time, is again qualitatively similar to that observed in the site 
approximation. 



z 



k 



(31) 
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V. DISCUSSION 



We studied the quasi-stationary properties of the DKCA in the one- and two-site approx- 
imations. Our study represents and extension of QS analysis, applied to continuous-time 
models exhibiting an absorbing-state phase transition in Ref . [^Of , to discrete-time processes. 



Compared with continuous-time processes, the numerical analysis of a discrete-time sys- 
tem is simpler, since it involves iteration rather than integration of a set of differential 
equations. While this is evident at the one-site level, at higher levels of approximation the 
advantage is tempered by the fact that starting from a given configuration, transitions to 
many (or all) other configurations are possible. The resulting need for combinatorial factors 
(e.g., T(N, Z; L)), complicates the analysis. 

An interesting result of our study is that the scaling behavior along the critical line is 
the same for the continuous-time contact process (and the closely related Malthus-Verhulst 
process) as for the discrete-time DKCA. In particular, the QS order parameter decreases 
~ 1/VL in both cases, while the QS lifetime grows ~ \TL. While the universality of 
global scaling could have been anticipated on the basis of the central limit theorem, the 
similarity extends further, to include the detailed form of the scaling function governing the 
QS probability distribution and its associated moments. Thus the situation is analogous 
to that found numerically in studies of absorbing-state phase transitions: not only critical 
exponents, but moment ratios of the order parameter take universal values at the critical 
point 
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TABLES 



TABLE I. DKCA transition probabilities 


<7iK'-i,0i +1 


1,1 1,0 


0,1 


0,0 


1 


P2 Pi 


Pi 








l-p 2 I- Pi 


I -Pi 


1 


TABLE II. Allowed values for Z on a ring. 


N 


z 






0,1 









2,..., L/2 


0,..., N-l 






L/2, ...L — 1 


2N - L,...,N - 1 






L 


L 
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Appendix 

To evaluate T(N, Z; L), the number of configurations on a ring of L sites with exactly iV 
active sites and Z nearest-neighbor pairs of active sites, we observe that the associated gen- 
erating function ((x, y; L) = J2n,z r(iV, Z; L)x z y N , can be written as the partition function 
for a one-dimensional lattice gas: 

C(x,y;L) = E ••• E x5>"*" y5>* , (35) 

0-1=0 cr L =0 

with o"at + i = o\. (We note in passing that x = e@ J and y = for the lattice gas with 
nearest-neighbor interaction J, chemical potential /x, and inverse temperature (3.) The 
partition function is evaluated using the transfer matrix T(er, a') = x aa ' y( a+a '^ 2 : 

((x,y;L) = TrT L (36) 

= Af + \ L 2 , (37) 

where X 1:2 are the eigenvalues of T: 

Ai, 2 = l/2 (l+xy±y/(l-x y y + 4y) . (38) 



For L even, we have: 

L/2 



{a + b) L + {a-b) L = 2Y J ^a 2n b 



2nvL—2n 



(39) 



so that 



leading to 



Af + A , . ^ g f. ^ (l\ M (in] H (z«! . (41) 



n=0 m=0 p=0 <jr=0 

The coefficient of x z y N is: 

L/2 2n / r \ / \ / L 



2n I \m I \ p I \ q I (Ay) 



(42) 

The above expression is evaluated numerically. 
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FIGURE CAPTIONS 



FIG. 1. DKCA phase diagram - simulation results from Ref. [21 



FIG. 2. Evolution of the probability distribution, conditioned on survival, in the one-site 
approximation. The initial distribution is po(N) = 6^,50- System size L = 100/ pi=p2 = 0.5. 



FIG. 3. Evolution toward the QS state in the one-site approximation; parameters as in Fig. 
2. Upper left: decay rate 7; right: mean number of active sites. Lower left: moment ratio 
m; right: skewness and kurtosis. 

FIG. 4. QS distribution at several different points in the DKCA phase diagram, for L = 100. 



FIG. 5. Quasi-stationary density for the DKCA, in the site approximation, for several system 
sizes. The inset shows moment ratio m versus pi for the same system sizes. 

FIG. 6. Moment ratio m versus L~ x l 2 at the critical point. □: one-site approximation, 
Pi=P2 = l/2;°: pair approximation, pi = 0.6306 , p 2 = 0.5. 

FIG. 7. Scaling plot of the QS probability distribution in the one-site approximation at the 
critical point {pi=P2 = 1/2), for L = 10 3 , 10 4 , and 10 5 (curves with maxima approaching 
the y-axis as L increases), and the asymptotic scaling function for the contact process found 
in Ref. J20| (with maximum at x = 0). Inset: Scaling plot of the QS distribution as in Fig. 
7, for L = 2000, p 1 = 1/2, and p 2 = 0.25, 0.5, 0.75, 0.9, and 0.999. The first four curves 
collapse, while the last has a broader distribution. 

FIG. 8. Quasi-stationary active-site density versus p\ for the DKCA in the pair approxi- 
mation, for several system sizes. The inset shows the moment ratio m for the same system 
sizes. 

FIG. 9. QS pair density probability distribution, conditioned on survival, at criticality 
(pa = 0.5, p x = 0.6306). 
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